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Two-particle correlations in relative rapidity and azimuth are studied for the p-Pb collisions at the LHC energy of 
y^SNN = 5.02 TeV in the framework of event-by-event 3 + 1-dimensional viscous hydrodynamics. It is found that for 
the highest-multiplicity events the observed ridge structures appear in a natural way, suggesting that collective flow 
may be an important element in the evolution of the system. We also discuss the role of the charge balancing and the 
transverse- momentum conservation. 
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Recently, the CMS collaboration has measured the two- 
dimensional (2D) two-particle correlations in relative pseu- 
dorapidity and azimuth for the proton-Pb (pPb) collisions 
at the energy of ^/snn — 5.02 TeV observing long- 
range structures, in particular, the near-side ridge. This 
fact contributes significantly to the on-going discussion on 
the very nature of the dynamics and its potential collec- 
tivity for high-multiplicity systems created in relativistic 
proton-nucleus and proton-proton (pp) collisions. We re- 
call that in relativistic nucleus-nucleus collisions the ridge 
structures are naturally explained by the presence of 
the collective flow arising from hydrodynamics. In [B[ one 
of us (PB) has studied the hydrodynamic evolution in the 
most central high-energy pPb and deuteron-Pb collisions, 
with the conclusion that a sizable elliptic and triangular 
flow can be formed. 

Some degree of collectivity has even been suggested for 
pp collisions [6W12I] . In high multiplicity pp events a same- 
side ridge is observed in the 2D correlations functions [l3[ • 
This feature could signal the presence of azimuthal corre- 
lations in the gluon emission in the initial state 
The same-side ridge observed in pp collisions could also 
result from the azimuthal asymmetry in the collective ex- 
pansion of the small droplet of matter created in the reac- 
tion [Hill. 

The hydrodynamic picture in the most central pPb colli- 
sions at the LHC is better justified that in pp interactions. 
The size of the system is comparable to the case of most pe- 
ripheral Pb-Pb collisions, where the elliptic flow has been 
observed. The initial density profile and its azimuthal 
asymmetry in pPb collisions can be estimated with the 
models tested in nucleus-nucleus reactions, whereas the 



shape of the small fireball that could be formed in pp col- 
lisions is less under theoretical control. In this Letter we 
model the 2D correlations for the most central pPb colli- 
sions in the CMS setup, reproducing the basic features of 
the data, in particular the presence of the near-side ridge. 
The hydrodynamic description of the fireball, resulting in 
collective flow, may therefore be an efficient approach even 
for small colliding systems, bringing in, among other pos- 
sible sources, a sizable component to the 2D correlations. 

The basic object of our study is the two-particle cor- 
relation function in relative pseudorapidity and azimuth, 
defined as |l| 
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where Arj and A(j) are the relative pseudorapidity and az- 
imuth of the particles in the pair. The signal is defined 
with the pairs from the same event, 
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while the mixed-event background distribution is 
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The number of particles N (denoted by CMS as iV tr i g ), 
is the number of charged particles in a given centrality 
class and acceptance bin, corrected for the experimental 
efficiency. The introduction of the central bin content, 
B(0, 0), brings in the interpretation of Eq. (p} as the av- 
erage number of correlated pairs per trigger particle. In 
our simulations we directly compute the right-hand side of 
Eq. ©. 

In pPb collisions, as the small interaction region fluc- 
tuates widely from event to event, one has to run the 
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(a) N > 110, 1<p <3GeV 




(a) N k >110, 0.1<p T <1GeV 



(b) 90<N, rk < 110, 1<p T <3GeV 
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Figure 1: The per-trigger-particle correlation function C tr i g (A77, A<j>) 
of Eq. |[T} for two most central ccntrality classes of the CMS Collabo- 
ration. The transverse momentum of each particle of the pair satisfies 
1.0 < p T < 3.0 GeV. 




Figure 2: The per-trigger-particle correlation function Ct r i g (Ar7, A</>) 
of Eq. JTJ for most central events, 0.1 < Pt < 1-0 GeV. The results 
including local charge conservation effects at the end of the hydrody- 
namic evolution are shown in panel (a), whereas panel (b) presents 
the case with non-flow correlations from the resonance decays only. 
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costly event-by-event simulations (e-by-e) @, Il9rl26j |. 
It is well known that the inclusion of the e-by-e fluctu- 
ations is important for the proper description of the initial 
eccentricity and triangularly, translating into the elliptic 
and triangular flow [27H30l | . For the small pPb system 
viscosity plays an important role |3lj . moreover, the den- 
sities are strongly rapidity dependent, hence viscous 3 + 1- 
dimensional (la . [32[ hydrodynamics must necessarily be 
used. 

We use the hydrodynamic model as described in [B[ to 
model the pPb system at highest centralities. The ini- 
tial condition is generated with GLISSANDO [33T| , imple - 
menting various variants of the Glauber model [34437 1 . 
The parameters of the calculations are similar as in [5|, 
except that they are adjusted for the collisions energy of 
y^SNN = 5.02 TeV. The nucleon-nucleon cross section is 
67.7 mb, moreover, we use a realistic (Gaussian) wound- 
ing profile [38[ for the NN collisions. To provide suitable 
initial conditions for the hydrodynamic evolution, the po- 
sitions of the participants are smoothed with a Gaussian 
with the width of 0.4 fm. 

We use the following initial profile in the space-time ra- 
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with 7/0 = 2.5 and a n — 1.4. The starting time of hydrody- 
namics is r = 0.6 fm, and the ratio of the shear viscosity to 
entropy density is rj/s — 0.08. The expected multiplicity in 
central pPb collisions is extrapolated linearly in the num- 
ber of participant nucleons from the minimum bias results 
of the ALICE collaboration [3{|. Accordingly, the initial 
entropy per participant in the fireball is adjusted. For ev- 
ery centrality we produce 300 initial configurations that 
are evolved with hydrodynamics to obtain freeze-out hy- 
persurfaces of constant temperature Tf = 150 MeV. Then, 
for each freeze-out configuration we generate 1000 THER- 
MINATOR events to efficiently improve the statistics. 

In the statistical emission model the non-flow correla- 
tions from resonance decays are built in. Additional cor- 
relations can appear due to local charge conservation [40( . 
Observation indicate that this charge balancing happens 
at hadronization 4T43j], i.e. at the late stage of the evo- 
lution. The mechanism creates a significant contribution 
to the 2D correlations functions [4| and is included in the 
simulations presented below. Pairs of opposite-charge par- 
ticles and their antiparticles are emitted locally at freeze- 
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out, with a thermal spread in their relative momenta. 
Another important source of correlations comes from the 
global transverse-momentum conservation 44, 45] . We im- 
pose approximately this constraint by requiring that the 
sum over the particles in the generated event fulfills the 
condition 
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We have found numerically that limiting the total trans- 
verse momentum to Pt — 5 GeV is sufficient; further 
reduction does not affect the studied quantities. This 
amounts for retaining about 8% of the least-Pr events from 
our sample. 

We apply the hydrodynamic model to the two most cen- 
tral centrality classes used by the CMS Collaboration. The 
centrality of the events is defined based on the charged 
particle multiplicity in the CMS acceptance. A good ap- 
proximation of the centrality cuts in our model is repre- 
sented by simple conditions on the number of the partici- 
pant nucleons. The most central collisions with N part > 18 
amount to 3.4% of most central events in the Glauber 
Monte Carlo model. The second most central class is de- 
fined by 16 < Np ar t < 17 and sums up 4.4% of the cross 
section. Cuts on the final multiplicity in the calculations 
instead of N part could be used, once the model of the ini- 
tial state were supplemented with effects of fluctuations of 
the energy deposited in each elementary collision 46-49j|. 
but this is not crucial for our study. 

In the hydrodynamic model the multiplicity fluctuations 
are largely decoupled from the collective expansion phase. 
Our model gives realistic predictions on the collective flow, 
but the multiplicity distribution cannot be reliably calcu- 
lated. This has a consequence for the normalization of 
the correlation functions. By integrating the per trigger 
correlation function one obtains 



i.e., the ratio of the average number of pairs over average 
multiplicity in a given acceptance window. In the presence 
of correlations from collective flow only, a more robust 
observable is the 2D correlation function normalized by 
the number of pairs instead of N in Eq. JTJ) . 

In the hydrodynamic model collective flow dominates 
in the correlation function for 2 < \Ar/\ < 4. There- 
fore to make a meaningful estimate of the hydrodynamic 
component in the 2D correlation function, we rescale the 
calculated functions to get the same subtraction constant 
Czyam in the zero-yield-at-minimum (ZYAM) procedure. 
We use the ZYAM values as quoted by the CMS Collabo- 
ration for each multiplicity and px bin [1]. Such rescaled 
correlation functions, called normalized correlation func- 
tions in the following, should be used to estimate the con- 
tribution of the collective flow to the ridge observed in the 
experiment. 
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0.1 <p T <1.0 GeV 
90<N, rk <110 



1.0<p T <2.0 GeV 
• CMS pPb 5.02TeV 




1 , 2 3 

|A *| 



1 |A*| 2 3 



Figure 3: The projected and ZYAM-subtracted correlation function 
in the region 2 < |Ar/| < 4 for the two most central bins in multiplic- 
ity (panels extending horizontally) and two pt intervals (panels ex- 
tending vertically) for the pPb collisions. The CMS measurement [l[ 
is shown as dots. The results of our hydrodynamic model with the 
normalized correlation functions are shown with the solid lines. The 
dashed lines show the results of the hydrodynamic model with sub- 
traction of the model ZYAM values and no rescaling. 



In the following we describe the results obtained with 
our simulations. We begin with the correlation function of 
Eq. shown in Fig. [T] for the most central collisions with 
two different cuts imposed on the transverse momentum 
of each particle in the pair. The 2D correlations func- 
tion presents similar features as the experimental one [![. 
A sharp same-side peak is formed due to the resonance 
decays and the local charge conservation QJ. The ob- 
served additional correlations from jet fragmentation at 
small A<fi-Ari, or the Bose-Einstein and Coulomb corre- 
lations, are not included in our model. The away- and 
same-side ridges are formed in the whole range of Aij. 
The shape of these ridges is determined mainly by the first 
3 harmonics in the relative azimuth. The first harmonic 
comes predominantly from the transverse-momentum con- 
servation and is seen as a tendency for the back-to-back 
emission. The second and third harmonics are provided 
by the collective expansion of the initial fluctuating source 
and describe well the shape and the width of the same- and 
away-side ridges. As expected 11, IH, the- collective ellip- 
tic flow leads to the formation of the same-side ridge in the 
2D correlation functions, which is our basic observation. 

A qualitatively different behavior is visualized in Fig. [2l 
At low pt the correlation displays a ridge (panel a) in the 
azimuthal angle direction (near Ai] — 0), which is due to 
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charge balancing in the hadronization and resonance de- 
cays. Since the harmonic flow is small at low px, only 
weak traces of the ridges (near Atfr — and A<f> — it) are 
visible. To evidence the effect of the local charge conser- 
vation in the formation of the A77 ~ ridge, we show in 
panel b the correlation obtained without charge balancing. 
The structure is much less pronounced now, as it is due to 
resonance decays only. We note, interestingly, that simi- 
lar structures, with a ridge in the A(f> direction in the 2D 
correlation functions, have been observed in pp collisions 
at 7 TeV by the ALICE Collaboration 

In the higher px bins (Fig. [TJ, the average transverse 
flow forms prominent ridges at A<f> ~ and A(f> ~ n which 
hide the A 77 ~ ridge. Also, as the flow increases, the 
unlike -sig n pairs become collimated in the azimuthal an- 
gle 0, l5ll | which contributes to the rise of the central peak. 

In a small system, a substantial part of particles is ex- 
pected to be emitted from the corona, without subsequent 
rescattering [52H55]. Particles emitted from the corona 
give a separate contribution to the 2D correlation function. 
In particular, one expects less collimation from charge bal- 
ancing [4(| • The shape of the Ary ~ ridge could be used 
to separate the contributions from the thermalized core 
and the corona in pPb reactions. In the kinematic re- 
gion IA77I > 2, selected for the analysis of the projected 
correlation function Ct r i g (A77, A<f>), the short-range charge 
balancing effects are not important. 

Next, in order to compare quantitatively to the CMS 
data we show in Fig. [3] the averaged correlation func- 
tion 

1 rfiV pair r r 

dArjC^A^A^)/ dA V . (7) 

iV dlXtp J2<|Ar,|<4 J2<|Ar;|<4 

We note that our e-by-e hydrodynamic simulations (lines) 
have the desired two-ridge structure, which is generated 
by the harmonic flow. The incorporated transverse mo- 
mentum conservation increases the relative strength of the 
away-side ridge. The normalization of the correlation func- 
tion is chosen to reproduce the normalization in the ex- 
periment. The function Ct r i g (A?7, A<f>) is thus rescaled to 
obtain the same value of the parameters Czyam as for the 
data points in Fig. 2 of The normalization procedure 
assures that the ratio (N(N—1)) > / < (N > is the similar 
as in the experiment. These normalized results, denoted 
with the solid lines in Fig. [3] are in good agreement for the 
1.0 < pr < 2.0 GeV case, and reproduce part of the ridge 
amplitudes in the lowest pr bin, 0.1 < pr < 1-0 GeV. The 
deviations may appear for many different reasons, namely, 
other sources of correlations are present, the initial ec- 
centricities calculated in the Glauber model without fluc- 
tuation at sub-nuclear scales may be underestimated, or 
contributions from the non-thermal corona in the inter- 
action region are important. In Fig. [3] are also shown 
the results obtained without normalizing the correlations 
functions but just subtracting the values at the minimum 
of the yield (dashed lines), Czyam; in that case the pa- 
rameters Czyam are a factor 1.1 to 2.5 larger than in the 
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Figure 4: The flow coefficients ^(Ajj) (upper lines) and v:j(Ari) 
(lower lines) calculated from the two-particle correlations as function 
of the relative pseudorapidity of the particles in the pair. The solid 
and dashed lines are for the unlike- and like-sign pairs, respectively. 
The central peak is due to charge balancing and, to a lesser extent, 
resonance decays. 

experiment. Thus the better agreement with the experi- 
mental points of the dashed lines is partly accidental, due 
to a mismatch in the particle multiplicities. Nevertheless, 
we show these curves, as the issues related to proper nor- 
malization between experiment and a model are far from 
trivial. 

In view of the recent results of the CMS Collaboration 
for the 2D correlations functions in pPb collisions, it is 
interesting to look at the possibility of measuring directly 
the harmonic flow coefficients. We plot the elliptic and tri- 
angular flow coefficients as function of the pseudorapidity 
gap in Fig. [4] The quantities are obtained in our hydrody- 
namic model from the Fourier decomposition of the cor- 
relation function Ct r i g (A?7, A(f>) [4]. The non-flow effects 
present in our model are important only for pairs of small 
pseudorapidity separation. In the intervals 2 < \Ar/\ < 4 
the non-flow effects from the resonance decays and the lo- 
cal charge conservation can be neglected. We note that 
the flow coefficients in Fig. [4] are sizable, thus could be 
measured. It must be noted, however, that other sources 
of non-flow correlations may be present also in that kine- 
matic region, but with smaller amplitudes, as measured in 
pp collisions [IS ]. 

In this Letter we have explore the possibility that the 
azimuthally asymmetric collective flow is generated in the 
hydrodynamic expansion of the fireball created in pPb 
collisions 0. The flow asymmetry together with the 
transverse-momentum correlations is capable of reproduc- 
ing the observed form of the 2D correlation functions in 
A-q-Acj) for the studied two highest multiplicity and two 
lowest pt bins, where hydrodynamics is expected to ap- 
ply. The agreement is semiquantitative, but very sugges- 
tive. The found important contribution of correlations of 
hydrodynamic origin (flow) does not exclude the presence 
of other sources of correlations between emitted particles, 
such as jets, string decays, or correlated gluon emission. 
We also note that the hydrodynamic description can be im- 
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proved in many ways, including fluctuations at sub-nuclear 
scales, the early flow, contributions from the corona, or 
varying the viscosity coefficients. The significant increase 
of the ridge amplitude when going from pp to pPb colli- 
sions indicates that the harmonic coefficients of the collec- 
tive flow become sizable and could be directly measured 
in spite of the background of the non-flow correlations in 
the small system. 

We notice the formation of the ridge at Ai] ~ at low 
transverse momenta in our pPb simulations, similarly to 
the recent findings of [50]. 
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